arXiv: 1504.00697v2 [quant-ph] 15 Jan 2016 


Classically entangled optical beams for high-speed kinematic sensing 

Stefan Berg-Johansen,^’^ Falk Toppel,^’^’^ Birgit Stiller,^’^ Peter Banzer,^’^’^ Marco Ornigotti,^ 

Elisabeth Giacobino/’^ Gerd Leuchs/’^’^ Andrea Aiello/’^ and Christoph Marquardt^’^ 

^ Max Planck Institute for the Science of Light, Gunther-Scharowsky-Strafle 1/Bldg. 24, 91058 Erlangen, Germany 
^ Institute for Optics, Information and Photonics, Universitdt Erlangen-Niirnherg, Staudtstrafie 7/B2, 91058 Erlangen, Germany 
^Erlangen Graduate School in Advanced Optical Technologies (SAOT), Paul-Gordan-Strafie 6, 91052 Erlangen, Germany 
Department of Physics, University of Ottawa, 25 Templeton, Ottawa, Ontario, KIN 6N5 Ganada 
^Institute of Applied Physics, Priedrich-Schiller University Jena, Max-Wien Platz 1, 07743 Jena, Germany 
^ Laboratoire Kastler Brossel, Universite Pierre et Marie Gurie, 

Ecole Normals Superieure, GNRS, 4 place Jussieu, 75252 Paris, Prance 
(Dated: January 18, 2016) 


Tracking the kinematics of fast-moving objects 
is an important diagnostic tool for science and 
engineering. Existing optical methods include 
high-speed CCD/CMOS imaging [T], streak cam¬ 
eras [2], lidar [3], serial time-encoded imaging [4] 
and sequentially timed all-optical mapping 
Here, we demonstrate an entirely new approach 
to positional and directional sensing based on the 
concept of classical entanglement [BHS] in vector 
beams of light. The measurement principle re¬ 
lies on the intrinsic correlations existing in such 
beams between transverse spatial modes and po¬ 
larization. The latter can be determined from 
intensity measurements with only a few fast pho¬ 
todiodes, greatly outperforming the bandwidth 
of current CCD/CMOS devices. In this way, our 
setup enables two-dimensional real-time sensing 
with temporal resolution in the GHz range. We 
expect the concept to open up new directions in 
metrology and sensing. 

Vector beams of light with cylindrical, non-uniform po¬ 
larization patterns [9] have found application in diverse 
areas of optics such as improved focusing m, laser ma¬ 
chining m, plasmon excitation [12] , metrology [13] , op¬ 
tical trapping m and nano-optics [Mill- Recently, vec¬ 
tor beams have attracted attention [T8H22] due to a sim¬ 
ple but striking property: when viewed as a superposition 
of transverse electromagnetic modes with orthogonal lin¬ 
ear polarizations, the nonseparable mode function of a 
radially polarized vector beam is mathematically equiv¬ 
alent to a maximally entangled Bell state of two qubits 
known from quantum mechanics. In contrast with the 
canonical Bell states in quantum optics, where two pho¬ 
tons are entangled in polarization and exhibit non-local 
correlations when spatially separated, this “classical en¬ 
tanglement” in vector beams is necessarily local as it ex¬ 
ists only between different degrees of freedom of one and 
the same physical system [23] . 

However, these correlations have recently been shown 
to represent a valuable resource. Vector beams have 
been used to violate an analogue of Bell’s inequality for 
spin-orbit modes na HO] and have led to continuous- 
variable entanglement between different degrees of free¬ 
dom [24]. In addition, vector beams have been used 



FIG. 1. Classical entanglement, a. (top) The trans¬ 
verse electric held distribution of a radially polarized beam 
of light. The orange density plot shows the beam’s doughnut¬ 
shaped intensity distribution, while black arrows indicate the 
position-dependent instantaneous local direction of the elec¬ 
tric held vector, (bottom) The beam’s global polarization 
state is shown in a Poincare sphere representation. Initially, 
the light held is globally unpolarized, with the Poincare vector 
located at the origin, b. If an opaque obstacle is brought into 
the beam, the global polarization takes on non-zero values 
according to the obstacle’s position within the beam. This 
method allows the object’s kinematics to be inferred from a 
polarization measurement alone. 


to implement classical counterparts of quantum proto¬ 
cols [25] [26] . Promising proposals include an application 
to the study of quantum random walks m and real¬ 
time single-shot Mueller matrix measurements [28], and 
a scheme for measuring the depolarization strength of 
a material has been implemented [29]. In the present 
work, we demonstrate for the hrst time a fully opera¬ 
tional application of classical entanglement to high-speed 
kinematic sensing. 

Several techniques are nowadays available for sens¬ 
ing the kinematics of fast-moving objects m- Each 
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comes with its own strengths and drawbacks. For exam¬ 
ple, high-speed imaging is typically limited to capturing 
only a small number of frames, while pump-probe tech¬ 
niques require the recorded event to be repeated identi¬ 
cally many times. Ideally, one seeks a solution that is 
capable of performing fast sensing continuously, in real¬ 
time and from a simple setup, employing only standard 
equipment and offering flexibility in the choice of wave¬ 
length. By using the nonseparable mode structure of 
cylindrically polarized beams (see Figure [^, one only 
needs to detect changes in polarization, thus fulfilling all 
the above requirements at the same time. In the follow¬ 
ing, we first discuss the physics of vector beams and then 
introduce the technique of sensing and show the results 
of our experimental investigations. 

The electric field of a general non-uniformly polarized 
paraxial beam can be written as: 


E{p,z) =eifi{p,z) + e 2 f 2 {p,z), (1) 


where the vectors ei, €2 determine the beam polariza¬ 
tion, the scalar functions /i(p, 2 :), /2(p, 2 :) set the wave- 
front and p = XX ^ yy is the transverse position vec¬ 
tor (see “Methods” below). The expression 0 is non- 
separahle^ namely it is not possible to rewrite it as the 
simple product of only one polarization vector and a 
single scalar function. In this sense, Eq. 0 has the 
same mathematical structure as a two-qubit entangled 
quantum state [23]. It is a well-established result of 
mathematical physics that any two-dimensional field of 
the form Q can be recast in the so-called Schmidt 
form F^(p,T) = \/Ai UiVi (p, z) -h ^/X2 U2V2 (p, z) where 
{'Ui,'U 2 } and {vi^V 2 } form complete orthonormal bases 
in the polarization and spatial mode vector spaces, re¬ 
spectively, with Ai > A 2 > 0. If either Ai = 0 or A 2 = 0, 
the expression of E{p^z) is factorable and the beam is 
uniformly polarized. Vice versa, if A 1 A 2 7 ^ 0, the beam 
displays a non-uniform polarization pattern and is said 
to be “classically entangled”. Thus, analogously to a 
bona fide quantum state, in a non-uniformly polarized 
beam, polarization and spatial degrees of freedom are so 
strongly correlated that if, by any means, one alters the 
beam’s spatial profile, then the polarization changes ac¬ 
cordingly. Our sensing technique relies precisely upon 
this peculiar phenomenon. Owing to the classical entan¬ 
glement exhibited by the beam, we are able to retrieve 
information about the position of a moving object par¬ 
tially obstructing the beam only by measuring the polar¬ 
ization of the latter: no spatially resolving measurements 
are needed. Since polarization measurements can be per¬ 
formed at GHz rates, with our system we are able to track 
very fast objects. 

For a field -E(p, z) in the Schmidt form, the measurable 


Stokes parameters can be written as 


« 5 o — Ai + A 2 , ( 2 a) 

51 = (Ai - A 2 ) - \ayf), (2b) 

52 = (Ai - A 2 ) (a^a* -h alay ), (2c) 

53 ~ '^(Ai A2) {cLx^y ^x^y) f 

where ax = ui - x and ay = Ui - y. In a radially polarized 

beam one has Ai = A 2 and thus 5 i = 52 = 53 = 0 , 


reflecting the fact that such a beam appears completely 
unpolarized in the absence of an obstruction. 

When an opaque object cuts across a non-uniformly 
polarized beam, the spatial and polarization patterns of 
the latter vary with time according to the obstructing ob¬ 
ject’s instantaneous position, as described by its central 
coordinates xo{t),yo{t). It is not difficult to show that 
for such a modified beam, Eqs. 0 are still valid, pro¬ 
vided that Ai, A2, ax, ay are regarded now as functions of 
xo{t)^yo{t)' When the values of the Stokes parameters 
5 o, 5 i, 82,53 are replaced by the measured ones on the 
left side of Eqs. 0 , these can be regarded as a nonlinear 
algebraic system of four equations for the two variables 
xo(t),^o( 0 : which can be solved by means of suitable al¬ 
gorithms. In this way, the instantaneous trajectory of 
the object is recovered. 

The experimental setup is shown in Eigurej^ We pre¬ 
pare a continuous-wave laser beam in a radially polarized 
mode. The beam impinges on a moving sample. Sub¬ 
sequently, half-waveplates and polarizing beam splitters 
are used to project the beam onto its horizontal, verti¬ 
cal, diagonal and anti-diagonal polarization components. 
Einally, a network of four InGaAs photodetectors with 
4 GHz 3 dB-bandwidth measures the individual projec¬ 
tions, from which the Stokes parameters 80 , 81 and 82 
can be straightforwardly obtained (see “Methods”). Eor 
the particular case of a radially polarized mode, the 83 
parameter is always zero. An auxiliary camera is used for 
additional visual verification and beam characterization. 

In order to demonstrate the system’s broad applicabil¬ 
ity, three types of measurement are carried out. Eirst, a 
metal rotor is made to turn about the beam axis (Eig- 
ure § 1 ). By sampling the Stokes parameters during the 
motion of the rotor, the instantaneous value of its an¬ 
gle of rotation Oq is succesfully inferred. An accuracy of 
4.1° (mean error) is achieved without correcting for beam 
imperfections and detector coupling. 

Second, a metal sphere is moved across the beam (Eig- 
ure 0 )- We measure the Stokes parameters with an ac¬ 
quisition time of 250 ps at each position. A Bayesian algo¬ 
rithm is used to estimate the sphere’s position from these 
data (see Supplementary Material). The inferred trajec¬ 
tory shown in Eigurej^ is seen to be in good agreement 
with the actual trajectory. As one expects, the inference 
is particularly successful in areas where the beam has a 
high intensity, i.e. where the Stokes parameter modula¬ 
tion introduced by the sphere results in a higher signal- 
to-noise ratio. 
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FIG. 2. Experimental setup, a. A continuous-wave laser beam is prepared in a radially polarized mode by a liquid crystal 
mode converter, b. The beam impinges on an opaque object whose motion in space modulates the beam’s polarization Stokes 
parameters, c. A polarization-independent beamsplitter (BSl) taps off 10 % of the beam for inspection by a conventional camera. 
This allows for mode pattern characterization and independent verification measurements, d. A polarization-independent 
50/50 beamsplitter (BS2) divides the beam up for projection onto its linear polarization components via a pair of polarizing 
beamsplitters (PBS) and a half-wave plate (A/2). The projections are simultaneously measured by four InGaAs detectors with 
4 GHz bandwidth. By linear combination of the projection signals, the beam’s Stokes parameters are obtained. Knowledge of 
the Stokes parameters allows the object’s instantaneous trajectory to be reconstructed (see Figure |^. 


Third and finally, the setup’s real-time capability is 
demonstrated by focusing the beam and measuring the 
Stokes parameters during the transit of a knife edge mov¬ 
ing at (27 ± 2) ms“^ across the focal plane. (The beam 
is sufficiently gently focused that it is not dominated by 
longitudinal field components at the waist, see “Meth¬ 
ods”.) As seen from the captured data in Figure]^ the 
transit takes only 92 ns, after which the beam is fully cov¬ 
ered. From the shape of the recorded traces, the knife 
edge’s direction of motion, horizontal in this case, can 
be inferred (up to a 180° rotation, see Supplementary 
Material). The event is captured as a sequence of single¬ 
shot measurements, requiring only a single occurrence. 
Furthermore, since the measurement is triggered on a 
change in sq, the particular instant of occurrence does 
not have to be known in advance. As the Stokes param¬ 
eters are captured continuously there is no dead time in 
this measurement. This result clearly demonstrates the 
measurement technique’s potential for high-speed kine¬ 
matic sensing. The technique allows for the use of bucket 
detectors rather than spatial detectors, and the measure¬ 
ment can be as fast as the detectors. With the analog 
bandwidth of 4 GHz available in our setup, we would be 
capable of resolving even sub-nanosecond motions. 

All three measurements confirm the setup’s ability to 
perform quantitatively meaningful kinematic sensing at 
very high temporal resolutions. We note that the mea¬ 
surement precision is subject to random error from the 
electronic detector noise at high bandwidths. This be¬ 
comes dominant in the regime where the measured sam¬ 


ple has only a small overlap with the beam (as seen in Fig¬ 
ure or when the sample covers the beam completely. 
Some applications, such as precision sensing of objects 
moving within a confined region, may benefit from using 
a beam with a nonzero 53 Stokes parameter. Such beams 
have been suggested for the investigation of small parti¬ 
cle scattering [30]. Although they require an additional 
photodiode pair, such beams avoid the zero of intensity 
at the origin. We note, however, that the classical en¬ 
tanglement of such a beam is not maximal, and that the 
correlations between polarization and position are there¬ 
fore necessarily weaker. 

In summary, we have demonstrated that the classical 
entanglement manifested by vector beams of light may 
be used to detect the kinematics of very fast objects 
with GHz temporal bandwidth. The method is possible 
because for cylindrically polarized beams, spatial infor¬ 
mation may be acquired by measurements on the polar¬ 
ization degree of freedom, owing to the classically en¬ 
tangled mode structure. The method presented requires 
only standard optical components which are commer¬ 
cially available at a wide range of optical wavelengths 
and can easily be extended to the microwave regime. 
It allows for continuous, real-time measurement of two- 
dimensional spatial information with unprecedented tem¬ 
poral resolution. We suggest that due to its simplicity, 
the method may even be employed in noisy environments 
such as free-space channels. For example, existing lidar 
technologies based on time-of-fiight measurements may 
be enhanced by the new method. On the microscale, fo- 








cused classically entangled modes may provide a new ap¬ 
proach to precision measurements, for example of Brown¬ 
ian motion in the ballistic regime in two dimensions [31] . 
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FIG. 3. Rotation sensing and position tracking, a. A metal rotor (width m = (0.79 zb 0.01) mm) turns about the centre 
of a radially polarized beam (width wi = (1.95 zb 0.10) mm). Due to the beam’s classically entangled structure, the rotation in 
space induces a sinusoidal oscillation of the beam’s Stokes parameters. Measurements of the so, si and S 2 Stokes parameters 
allow the instantaneous angle of rotation to be inferred. Each data point was obtained by integrating over 200 ns, so that 
electronic noise is averaged out to within the data point width. Dotted curves show the theoretically expected values under 
the assumption of an ideal mode function, b. (left) A metal sphere (diameter d = (1.00 zb 0.01) mm) traverses a radially 
polarized light beam (width W 2 = (2.84 zb 0.10) mm), (center) The Stokes parameters so, si and S 2 vary as a function of the 
sphere’s position. Solid lines show the expected Stokes parameters as obtained from simulation, (right) The sphere’s trajectory 
is inferred from the measured Stokes parameters. The sphere is moved gradually in discrete steps of 50 pm, providing a 
calibrated, reproducible reference motion. To allow for a realistic comparison with a fast object, the acquisition time at each 
point is only 250 ps. The blue contours show the combined Bayesian 68% credible region R, while the gray shadow shows the 
sphere’s dimensions to scale. The theoretical model used to obtain the theory and simulation curves for a and b, respectively, 
is detailed in the “Methods” section. The position tracking algorithm is described in detail in the Supplementary Material. In 
both plots. So is normalised to its initial value, while si and S 2 are normalised to the instantaneous value of sq. 
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FIG. 4. Real-time sensing. A metal knife edge of thickness (3 zb 2) jim cuts across a focused radially polarized mode 
(theoretically estimated width — (2.0 zb 0.5) pm) at (27 zb 2) ms“^. The plot shows a sequence of single-shot measurements 
of the beam’s so, si and S 2 Stokes parameters during the knife edge’s passage until the beam is fully covered (normalised to 
the initial power), with a total duration of 92 ns. The sampling resolution reaches up to 100 ps. From the measured traces, the 
knife edge’s direction of motion can be easily inferred up to a 180° rotation. 
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Methods 

Experimental setup. We employed a continuous-wave diode 
laser (Sumitomo SLT5411) with wavelength A = 1550 nm, ampli¬ 
fied by an erbium-doped fiber amplifier (OPREL OFA20-1231S) to 
a power of 50 mW. The initially Gaussian beam profile was con¬ 
verted to a radially polarized beam using a commercially available 
liquid crystal device (Arcoptix S.A.) with custom anti-refiection 
coatings. The converted beam was Fourier filtered with a 75 pm 
pinhole placed between a pair of f = 50 mm plano-convex lenses 
and passed two tilted waveplates in order to correct residual 
circular polarization components. The rotor was a Cu-coated 
wire of width (0.79 ± 0.01) mm, placed in a rotation mount which 
was actuated by a stepper motor. The sphere had diameter 
(1.00 ± 0.01) mm and was made from steel. It was attached to 
an uncoated glass substrate of thickness 1 mm by a meniscus of 
epoxy resin which was small enough to be neglected optically. The 
substrate was mounted in a vertical position on a linear precision 
stage (PI M-410.DG, minimum incremental motion 0.1pm) 
and moved in steps of 50 pm. The plane of object interaction 
was located at 50 cm after the mode converter. The distances 
between object-camera and object-detectors were 10 cm and 80 cm, 
respectively. The beam had a width of (2.84 ± 0.10) mm in the 
plane of object interaction, and a Rayleigh range of approximately 
2 m. All beam widths given here and in the main text are 90-10 
knife-edge widths. The detectors used to measure the horizontal, 
vertical, diagonal and anti-diagonal polarization components 
were two Soliton ET-3500, an Alphalas UPD-35 and an Alphalas 
UPD-70, respectively. The detectors were temporally matched 
to within 1.3 ns, with the remaining offset being compensated by 
the oscilloscope up to the measurement bandwidth. For the rotor 
measurement only, a pair of BPDV2120R fiber-coupled balanced 
photoreceivers with 43 GHz 3 dB-bandwidth from the company 
u^t were used. Tailored phase plates were used to convert the the 
linearly polarized Hermite-Gauss projection modes into approxi¬ 
mate fundamental modes in order to fiber-couple them to these 
detectors. For the real-time measurement, the beam was focused 
using a pair of aspheric lenses with NA 0.65. While it is known that 
radially polarized beams display a strong longitudinal electric field 
component on axis when focused tightly ca, such components 
are not yet predominant in the moderate focusing regime found 
here. Gomparative measurements with an azimuthally polarized 
beam, for which no longitudinal field components appear during 
focusing, showed very similar overall behavior and duration. All 
signals were recorded with a digital storage oscilloscope (Agilent 
DSO81204B in half-channel mode, 8 GHz analog bandwidth, 
sinc-x-interpolation disabled). An InGaAs camera with resolution 
320 X 256 px (Xenics XS) was used for beam inspection and for 
generating spatially resolved beam polarization tomographies 
prior to the actual sensing measurement. The Bayesian tracking 
algorithm described in the Supplementary Material used these to¬ 
mographies in order to improve its estimate for the object position. 


Theoretical model. In order to obtain the theoretical pre¬ 
dictions and simulated values shown in Figure we developed 
the following simple model. Gonsider a radially polarized beam 
of light of angular frequency ca, propagating along the 2 ;-axis 
of a Gartesian reference frame and polarized in the xy-p\ane. 
In the paraxial approximation its electric field takes the form 
£(r, t) = Re lE(p, z) ex.p{—iujt)] with time-independent mode func¬ 
tion E(p,z) (Jones vector). In the ideal case of a perfect radially 
polarized beam the mode function is given by 


E(p,z) 


1 

71 


[* V’io(p, z) + y V’oi(p, z)]. 


(3) 


Therein 'ipmn{p,z), with m,n G {0,1,2,...}, is the Hermite- 
Gaussian solution of the paraxial wave equation of order N = m+n, 
while p = XX + yy denotes the transverse position vector. Please 
note that Eq. § is a particular case of the general expression in 

Eq. (Q. 

When an opaque obstacle is placed into the beam, covering a 


region A in the xy-plane, the electric field directly behind the ob¬ 
stacle is given by 

E'(p,z) = ll-JAip)mp,z), (4) 

where 


^a{p) 


f 1 if p e A, 
if p ^ A, 


(5) 


denotes the transmission function of the aperture complementary 
to the obstruction (Babinet’s principle). 

We are interested in the perturbed beam’s intensity after passing 
a linear polarizer making an angle (f G {0°, 90°, 45°, 135°} with the 
x-axis, as the Stokes parameters so, si, S 2 are easily calculated from 
these intensities. Representing Eqs. § and as two-component 
Jones vectors and describing the action of the linear polarizer by 
the projection matrix 


TD _ 

cos^ Lp sin (p cos p 

G/? — 

sin p cos p sin^ p 

lAiv) = j 

\P^E’{p,z)\^d'^p 


= Jll-XA{p)]\P^E(p,z)\^d^p (6) 

for the beam’s intensity behind the polarizer. With the help of 
Eq. the Stokes parameters, measured in the plane of the ob¬ 
structing object, take the form 


so (.4) = 7^(0°)+ 7^(90°), 

(7a) 

SI (.4) = 7^(0°)-7^ (90°), 

(7b) 

s2(.4) = 7^(45°)-7^(135°), 

(7c) 


and are functions of M, the region covered by the object within 
the beam. In simple cases, these equations may allow for 
straightforward algebraic inversion. To obtain the theoretical 
predictions shown in Figure (rotor), we have used the time- 
independent mode function of the ideal radially polarized beam 
from Eq. The rotation angle can then be obtained analytically 
as ^0 = arcsin(si/s 2 ). The simulated reference in Figure & 
(tracking) relies on experimentally measured distributions of 
\PipE{p, z)]"^ that have been substituted into Eq. This allows 
for the numerical evaluation of a likelihood function L which 
is used by the position tracking algorithm (see Supplementary 
Material). 
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SUPPLEMENTARY INFORMATION 


Position tracking algorithm. Here we explain the algorithm 
used to infer the transverse central coordinates (a;o,yo) of a spher¬ 
ical object from measurements of the polarization Stokes parame¬ 
ters s = (so, 51,52). Among the advantages of a numerical imple¬ 
mentation is that experimental beam imperfections and electronic 
detector noise can be accounted for in a straightforward manner. 
Eqs. 0 in the “Methods” showed how to compute the individ¬ 
ual Stokes parameters as a function of the transverse position of 
an object covering a known region A. Before the actual sensing 
measurement, we perform this computation for each possible ob¬ 
ject position in discrete steps using the experimentally determined 
spatially resolved beam polarization tomographies \P^E{p,z)\'^,ei- 
fectively yielding a look-up table (xo,yo) —5 as visualised by the 
surfaces in Supplementary Figure [^. By taking into account the 
normally distributed, experimentally determined statistical noise of 
each detector, the look-up table can be turned into a likelihood dis¬ 
tribution L (s|(a:o, yo))- During the real-time sensing measurement, 
given a measured Stokes vector s, the position can be inferred by 
finding the coordinates {xQ,yo) maximizing this likelihood distri¬ 
bution (see Supplementary Figure]^). 

The rotational symmetry of the radially polarized mode implies a 
twofold ambiguity in determining the object’s position, since rotat¬ 
ing any object by 180° about the beam’s center leaves the Stokes 
parameters invariant. This ambiguity can be lifted by making a 
number of additional assumptions, such as specifying the half-plane 
where the particle is located at time t = 0 (initial position), and 
that its velocity is small compared to the sampling frequency (con¬ 
tinuity of position), or that its acceleration stays within physically 
motivated bounds (continuity of velocity). Such assumptions can 
easily be accommodated within the framework of Bayesian proba¬ 
bility theory (BPT) [von Toussaint, U. Rev. Mod. Phys. 83 , 943- 
999 (2011)] where they take the role of prior information. BPT 
also assigns a formal role to other assumptions such as the size 
of the object to be tracked, and provides the means to estimate 
such parameters when they are not known a priori. However, in 
our laboratory setting we observe that the beam’s symmetry is al¬ 
ready broken by small beam imperfections (including a slight, but 
not complete, four-fold symmetry in intensity). Interestingly, we 
found that taking into account this overall asymmetry helps to de¬ 
termine the trajectory unambiguously with high certainty without 
such additional assumptions as described before. If necessary, the 
symmetry could also be broken in more controlled ways, e.g. by 
introducing a second beam, or by using a beam with a less sym¬ 
metric polarization structure [Beckley, A. M. et al. Opt. Express 
18 , 10777-10785 (2010)]. 

We emphasize that the algorithm used is non-recursive and does 
not require any repeated iterations on the same data. Rather, it 
relies mainly on data which can be pre-computed in the form of 
look-up tables to any desired accuracy in a computationally effi¬ 
cient way. Hence, the algorithm is highly parallelizable and thus 
compatible with the requirements of high-speed real-time signal 
processing. The algorithm’s performance is rather insensitive to 
the assumed shape of the object: assuming a square of the same 
area instead of a circle leads to identical inferred positions (up to 
the width of the 68% credibility region). Changing the assumed ob¬ 
ject size generally leads to a radial shift of the inferred trajectory, 
but leaves the main features of the trajectory intact. 

For the measurements shown here, we accounted for diffraction 
by assuming an effective sphere diameter at the detector plane 
which is 35% larger than the sphere’s actual physical dimensions. 
It is likely that the accuracy can be further increased by incorpo¬ 
rating a more detailed model of diffraction. However, the actual 
impact of diffraction and the appropriate model to incorporate it 
will depend on the particular application under consideration. 
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FIG. 5. Position tracking (experimental data), a. For a known beam and object geometry, a look-up table 
{xQ^yo) S{o,i, 2 } may be computed from Eqs. Q, giving the expected Stokes parameters as a function of the object’s 

transverse position within the beam. One such look-up table, obtained from experimental data, is visualized here as a set of 
three surfaces. The measurement of a particular Stokes vector may be regarded as a set of planes intersecting these surfaces. 
The intersection of the corresponding contours (shown on the bottom) indicates the position of the object, b. The upper plane 
shows the likelihood distribution L(s|(xo, yo)) for the measured Stokes vector s. The maximum of this distribution corresponds 
to the inferred position. Note that although the theory for an ideal mode predicts two solutions, this symmetry may be broken 
in practice due to slight irregularities in the mode pattern, as demonstrated by this particular example. 











